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Abstract 

We report on the ALPHA Collaboration's lattice B-physics programme based on A^f -2 0(a) improved Wilson 
fermions and HQET, including all NLO effects in the inverse heavy quark mass, as well as non-perturbative renormal- 
ization and matching, to fix the parameters of the effective theory. Our simulations in large physical volume cover 3 
lattice spacings a x (0.08 -0.05) fm and pion masses down to 190 MeV to control continuum and chiral extrapolations. 
We present the status of results for the b-quark mass and the B(s)-meson decay constants, /b and /b^. 



1. B-physics and lattice QCD 

Lattice simulations of QCD have established as a sound 
tool to compute strong interaction effects for accurate 
phenomenology in heavy flavour physics. For B-meson 
weak decays, which constrain the CKM Unitarity Tri- 
angle, lattice QCD results for the involved low-energy 
hadronic matrix elements in conjunction with experi- 
mental studies decisively contribute to stringent tests of 
the self-consistency of the Standard Model and comple- 
ment direct searches for New Physics. Since the sig- 
nificance of such precision tests in the beauty sector is 
predominantly limited by theoretical uncertainties, lat- 
tice computations with an overall accuracy of a few % 
are highly desirable. Let us highlight the "Vub-puzzle" 
illustrated in Fig. [T] which has drawn the community's 
attention in the recent past. IVubl can be determined from 
inclusive semi-leptonic processes B — > Xy^v, from ex- 
clusive semi-leptonic B — > ntv decays and from the lep- 
tonic one, B — > tv. In the latter two, the hadronic form 
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Figure 1: Observed tension among different \Vub\-cleterminations fj^; 
± 1 o" bands are shown. 

factor f+iq^) and the B-meson decay constant /b enter, 
respectively, so that lattice QCD input is required to ex- 
tract IVubl- As currently there is a ~ 3cr tension between 
its two exclusive (semi-leptonic and leptonic) determi- 
nation^ as well as an inconsistency with the estimate 
from inclusive decays, precision lattice QCD calcula- 
tions can contribute to resolve this tension l|4j-i6J. 



'At the ICHEP 2012 Conference, the Belle Collaboration reported 
a new result for .!^(B — ► tv) \2 3 1, obtained on basis of a new data set 
using a more sophisticated tagging of the B. Taken alone, this would 
yield a value for IVubl compatible with the exclusive semi-leptonic de- 
termination. However, this result has not yet been confirmed by other 
collaborations, and more data and a careful study of all systematics 
are still required before drawing final conclusions. 
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2. Non-perturbative HQET 

The particular challenge of B-physics on the lattice lies 
in the many disparate scales, ranging from its inverse 
extent (as IR cutoff) over the hierarchy of differently 
flavoured hadron masses up to the inverse of the lattice 
spacing a (as UV cutoff), to be treated simultaneously 
in the numerical simulations. Since lattice sizes that are 
computationally manageable today have am\, > 1, dis- 
cretization eff'ects get most severe for heavy quark sys- 
tems with b-quarks and escape brute force simulations. 
Our approach to lattice B-physics is therefore based on 
the Heavy Quark Effective Theory (HQET) for the b- 
quark fTl, ISl, which consists in a systematic expansion 
of its QCD action and heavy-light correlation functions 
in Agco/wib 1 around the static limit (m\, — > oo). 

The Lagrangian entering the heavy quark field's lat- 
tice action 5hqet - a'^Hx ■^hqet(jc) in HQET at NLO, 
i.e., including 0(l/mb) terms, reads: 

-%QET = lAh-^OfAh - ^^kin(5kin - Wspint^spin , (1) 

with i/Th satisfying P+iph - lAh. P+ - The pa- 

rameters (jjkin and Wspin are formally 0{l/m\,) and mul- 
tiply the dimension-5 operators <9kin - lAhD^iAh and 
<5spin - 'l^h'^^ 'Ah. representing interaction terms due to 
the motion and the spin of the heavy quark. Thus, S hqet 
has ©(Ag^p/m^) truncation errors, and lattice artifacts 
only scale as (aAqcd)" rather than (am\,)". Analogously, 
local composite fields are introduced in the effective lat- 
tice theory. For instance, the NLO HQET expansion of 
the zero-momentum projected time component of the 
heavy-light axial vector current can be written as 

A«r(.o) = zf^^o^i:x[ArW + c:^'A("(x)], 

A^'^'(x) = ^Ai(x)yor5</'hW, 

A<;'(x) = iAiWr5r/3(^'-^')'AhW, (2) 

being the spatial co variant derivative. The relation 
/ps V^PS = ( 1 Ao,r(0) I PS(p = 0) > to the pseudoscalar 
decay constant will be used to calculate /b,,, below. 

HQET treats the 0(l/mb) interactions terms in ([T]i as 
local space-time insertions in static correlations func- 
tions. For correlators of some multi-local fields O and 
up to l/wb -corrections to the operator itself (irrelevant 
when spectral quantities are considered), this means 

iO) = <(9>stat + a*Yj[ ^^niOOi^nix)),^, (3) 
+ <y,spin<(5(9.spin(ji:)).stat } , 

where ((9)stat is the expectation value in the static theory. 

Still, for lattice HQET applications to lead to precise 
and controlled results, two issues had to be solved. 



1. ) The exponential growth of the noise-to-signal ra- 
tio in static-light correlation functions with Euclidean 
time, caused by the linear divergence in the binding en- 
ergy E"'^' of the static-light system, which is particularly 
severe for the Eichten-Hill action [9]. This is overcome 
by so-called "HYP-smeared" LIOJ discretizations of the 
static quark action, improving the statistical precision of 
the correlators substantially |Tr,T2l. 

2. ) Operator mixing in the effective theory induces 
UV power divergences in the lattice spacing that must 
be subtracted non-perturbatively: The formal definition 
of lattice HQET and its composite fields in ([T]) and (|2]i 
involves the (a priori free) effective couplings 

= { ^Mbare , In Z"*^'^'^ , ' , Wkin , t^spin } ■ (4) 

Here, the energy shift mbare is an additive mass renor- 
malization. It absorbs the 1 /a-divergence of the static 
energy, E^^'-, and a 1 /a^-divergence at 0(l/mb). Hence, 
a phenomenologically relevant predictive power of lat- 
tice HQET is only guaranteed, once these HQET pa- 
rameters (i) — {(AJi) have been fixed non-perturbatively 
such that no uncancelled power divergences in 1/a, 
which would remain in perturbation theory lfT3l . can 
preclude to take the continuum limit. 

A solution to 2.) was developed in fTT] and relies 
upon a non-perturbative matching of HQET and QCD 
infinite volume. The implementation of this strategy by 
our collaboration has led to NLO HQET computations 
of the b-quark mass, B-meson spectroscopy and decay 
constants in the quenched approximation (A^f = 0) ifTSl - 
fTS), as well as in the more realistic two-flavour theory 
111 9^422 1 . of which we give an overview in the following. 

2.1. General strategy 

The computational strategy of our approach lfT4l[T6ll22l . 
in which matching and renormalization are performed 
simultaneously and non-perturbativel}|^ is sketched in 
Fig. |2] The matching part is performed in a small 
volume of L\ ^ 0.4 fm, where owing to am\, <K 1 
numerical simulations with a relativistic b-quark are 
feasible. The bare HQET parameters w,- of the La- 
grangian and the time component of the heavy-light ax- 
ial current are fixed by imposing matching conditions 
<D"QET(^^g) 1 cD)QCD(^ 0) = Hm«^o<l>Q™(z,fl). The 



-As soon as l/wib-corrections are included, matching must be 
done non-perturbatively in order not to spoil the asymptotic conver- 
gence of the series. Otherwise, the perturbative truncation error from 
the matching coefficient of the static term becomes much larger than 
the power con'ections ~ AgcD/Hib of HQET, as mj, — » t». 
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Figure 2: /dea of iMice HQET computations for B-pbysics phe- 
nomenology via a non-perturbative detennination of HQET param- 
eters from small-volume QCD simulations. Tbe step scaling method 
makes contact to physically large volumes Loo. The whole construc- 
tion is such that the continuum limit can be taken for all pieces. 



quark mass dependence of (non-perturbatively renor- 
malized) QCD is inherited by the HQET parameters w,. 
It enters through the dimensionless variable z = L\M, 
where M is the renormaHzation group invariant (RGI) 
mass II20I . Then a recursive finite-size scahng step 
Li ^ L/i- 2Li is used to reach larger volumes and lat- 
tice spacings a, by which connection with phenomenol- 
ogy in Loo ^ max(2 fm, A lm„) can be made. As a result 
of II22I . the HQET parameters Q, co - co{z,a), absorb- 
ing the logarithmic and power divergences of HQET, are 
now non-perturbatively known for renormalized QCD 
quark masses from the charm to beyond the beauty re- 
gion (parameterised by z e {4,6,7,9, 11, 13, 15, 18,21)) 
and for a's corresponding to the bare gauge couplings 
of the available two-flavour configuration ensembles in 
large volume used to compute B-physics observables. 



2.2. Large-volume computations and techniques 

Our large-volume gauge configuration ensembles are 
characterised by the plaquette gauge action and a sea 
of A^f = 2 mass-degenerate non-perturbatively 0(fl) im- 
proved dynamical Wilson quarks. To be able to extrap- 
olate to the physical pion mass, several pseudoscalar 
(sea) masses in the range (190 < mps < 440) MeV 
with L mps > 4 are considered such that finite-volume 
effects are expected to be negUgible. Moreover, the con- 
figurations cover 3 yS-values (5.2,5.3,5.5) with lattice 
spacings a e {0.08 fm, 0.07 fm, 0.05 fm) L23J to control 
the extrapolation to the continuum Umit. 



For the numerical simulations to generate these two- 
flavour QCD configurations, we employ M. Liischer's 
implementation of the Hybrid Monte Carlo (HMC) al- 
gorithm with domain decomposition 11241 and an adap- 
tion [25], which combines the deflated SAP solver f26l- 
l28l with mass preconditioning ||29l . chronological in- 
version ||30]| and multiple time scale integration 113111321 . 
Large trajectory length 13311341 and long runs are used 
to ensure that our ensembles are not biased by the crit- 
ical slowing down of the QCD simulations. For a care- 
ful and conservative error estimation, a binned jackknife 
procedure is applied (which is being cross-checked by 
the method advocated in ll34l ). All large-volume con- 
figuration ensembles have been produced and are shared 
within the Coordinated Lattice Simulations (CLS) eff'ort 
by several lattice QCD teams in Europe ll35l . 

Our determination of the B-meson spectrum and de- 
cay constants in two-flavour QCD is based on their 
HQET expansions in terms of the known HQET pa- 
rameters and associated HQET energies and matrix el- 
ements at the static and 1 / njb-order The latter are ex- 
tracted from static-light correlation functions evaluated 
on the available large-volume CLS ensembles by solv- 
ing the Generalised Eigenvalue Problem (GEVP) dis- 
cussed and applied in 1361 1371 ITTl ITSl . allowing for a 
better control of excited state contaminations of the cor- 
relators. More specifically, the GEVP analysis amounts 
to compute a(N x A^)-matrix of correlators with the de- 
sired static and 0(l/mb) insertions, where each entry of 
the matrix corresponds to a different Gaussian smearing 
level l38l of the light quark field in the B-meson in- 
terpolating quark bilinear In these computations, vari- 
ance reduction in the light quark sector is achieved by 
stochastic all-to-all propagators (with 8 noise sources 
and full time-dilution) 1391 1171 . while for the static 
quark propagators, we use two variants of the HYP- 
smeared static actions, HYPl/2 10111121 . akeady men- 
tioned in Sect. |2] Solving the GEVP numerically gives 
rise to new estimators for effective energies and hadron- 
to-vacuum matrix elements, which converge faster with 
Euclidean time separation than standard ratios, since a 
larger gap governs the excited state corrections. I.e., 
corrections to ground state energies and matrix ele- 
ments fall off in t and fo as oc exp{-(£'A'+i - Ei)t} and 
oc exp {-{En+i - El) fo) X exp {-(La -Ei)(t- to)], re- 
spectively, where to < t < 2to and labels the A^* 
excited state (and = 3 in practice). Our final esti- 
mates are then obtained as plateau averages over ranges 
conservatively chosen by varying tj^in for fixed fmax such 



that OCe"*^"*'"*'^') ~ o-sys 
ing our systematic errors. 

toETlinKlslEollITl. 
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3 (Tstat, thereby minimis- 
For more details we refer 



3 



3. Results 

In this section we summarise the status of resuhs of our 
Nf - 2 B-physics project, as it was reported at summer 
conferences in 2012, see also ||40', '42). A final account 
of our work will be given later, once the full statistics of 
all CLS ensembles has been analysed |41 1. 

Our HQET energies and matrix elements extracted 
from static-light correlators split into two sets, one for 
the B-meson sector with the valence quark masses set 
equal to the CLS sea quark mass values, and another one 
for the Bs-meson sector with the valence quark tuned to 
the physical strange quark Il23l . corresponding to a par- 
tially quenched setup. 




Figure 3: Joint chiral and continuum extrapolation of the heavy-light 
pseudoscalar meson mass ^inNLOHQETtothe ansatz ^ for fixed 
z. Recall tliat the z-dependence originates from the initial matching 
step to finite-volume QCD determining the HQET parameters. (Blue, 
red and green points refer to/? = 5.2, 5.3 and 5.5, respectively, while 
filled/open symbols belong to the HYPl/2 static actions.) 



3.1. The b-quark's mass 

To begin with, we apply the non-perturbative results on 
the HQET parameters from the matching step together 
with the HQET energies from the large-volume GEVP 
analysis of the CLS ensembles to calculate the b-quark 
mass. To this end one writes down the NLO HQET ex- 
pansion (i.e., to first order in l/»Jb) of the heavy-light 
pseudoscalar meson mass as 

'«B = '«bare + + ^^kin^""" + W,p,„£^P'" . (5) 

Beside the dependence on the light pseudoscalar (sea) 
mass mps and the lattice spacing a of the CLS configu- 
rations, we also have to account for the apparent heavy 
quark mass (z) dependence of the HQET parameters. In 



fact, this z-dependence can now be exploited to fix the 
HQET parameters w,- = W/(z, a) for once by imposing 
the condition 



mjiiz, m„, a) \ = ^ = 5279.5 MeV fV\ (6) 



in the continuum, which defines the physical value of 
the b-quark mass at NLO of HQET. 

For given z, the l.h.s. is obtained by evaluating (|5]l 
for each value of the lattice spacing and sea quark mass 
of the CLS ensembles, followed by a global fit of 
(simultaneously for two variants of the aforementioned 
HYP-smeared actions, HYPn, n = 1, 2) to the ansatz for 
a combined chiral (mps — > m„ = in^'^ - 135 MeV U]) 
and continuum (a — > 0) extrapolation. 



in^{z, mps, a, HYP«) = 
3f 



Biz) + Cmf, 



(7) 



with = f„^'"^ = 130.4 MeV [1| and f = 0.51(2) Eg 
the B*B7r-coupling in the static approximation for the 
b-quark. These extrapolations to the physical point are 
shown in Fig. [3] 
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Figure 4: z-dependence ofm^ in the continuum limit and graphical 
solution of^, which determines the physical b-quark mass Zb- (z = 
L[ M denotes the dimensionless RGI heavy quark mass.) 



The solution of (|6| for a = as the physical con- 
dition defining the b-quark mass yields the dimension- 
less RGI b-quark mass Zb = L\M\, and is illustrated in 
Fig. [4] Converting Li to physical units, achieved via 
setting the lattice scale through /k in [23 J . and translat- 
ing with 4-/3-loop RG running of the coupling/mass to 
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the conventional MS scheme, we obtain as our resuh for 
the b-quark's mass in HQET at 0(l/mb) in the Nf - 2 
theory presented at this year's summer conferences: 



Zb = 13.34(33)(13), 



m^^iWib) = 4.22(10)(4),GeV. 



(8) 



The first error covers all statistical and systematic errors, 
including those from the GEVP analysis, the various ex- 
trapolations and the scale setting, while the second un- 
certainty of about 1% stems from the quark mass renor- 
malization in QCD, entering the finite-volume match- 
ing step II20I I22I . and has to be added in quadrature. 
Since we find the difference of the NLO HQET result 
(|8]l and the corresponding number in the static approxi- 
mation (LO HQET) to be very small, we conclude that 
the truncation error of 0{A^^j^/m^) to (jsj in the HQET 
expansion is negligible compared to our overall error. 

Our result (jS} also compares very well with other re- 
cent determinations and the value quoted by the Particle 
Data Group, see the (incomplete) compilation in Tab. [T] 
Note that some determinations claim very small errors, 
although they are based on perturbation theory or lattice 
data with heavy quark masses in lattice units close to 1. 



mb'''(mb)/GeV 


remarks, method 


ref. 


4.347(48) 


lattice, A'f = 0, NLO HQET 


ca 


4.22(11) 


lattice, A'f = 2, NLO HQET 


eq. jsj 


4.29(14) 


lattice, A'f = 2, extrapolation 


mi 


4.164(23) 


lattice, Nf = 3, extrapolation 




4.163(16) 


perturbation theory & data 




4.236(69) 


perturbation theory & QCD inputs 


da 


4.18(3) 


PDG average 2012 


m 



Table 1 : Compilation of some recent determinations of nib ■ tor the 
lattice results 1 44|/ and I45t . the former uses extrapolations of relativis- 
tic data around the charm to known static limits, while the latter em- 
ploys moments of current-current correlators with Highly Improved 
Staggered Quarks (HISQ) i481 extrapolated to the b-scale. 146^ and 
l47\l rely on QCD sum rules. For more details, see the cited references. 



After the determination of the value of the physical b- 
quark mass (and thus Zh) from m^, we can fix the HQET 
parameters to <y,(fl) = a)i{zb,a) and use those in any 
successive HQET computation of B-physics obervables. 



3.2. B-meson decay constants 

To determine the B-meson decay constants /b and /b,, 
we now combine the HQET parameters with the ma- 



trix elements resulting from the GEVP analysis. Dis- 
tinguishing the heavy-light Bj-meson decay constants, 
where the light flavour can be either a valence {- sea) 
quark flavour r = {u, d) = 1 or a valence strange one, 
r = s, their NLO HQET expansions in terms of the 
HQET parameters w, read 



Zf^^'' (l + bf'am, 



(9) 



As explained before, the are understood to be taken 
at the physical b-quark mass, Zb, and the p^, X e 
{stat,kin, spin, A*^'*), denote the previously extracted 
GEVP plateau values of the associated effective matrix 
elements of Sect. 



2.2 



The improvement coefficient bj^ 
is known to 1-loop perturbation theor}|^from |49|. 

Due to non-perturbative 0(fl) improvement employed 
in our computations, /b_ (with r = 1, s and /b, = /b), 
approaches the continuum limit quadratically in the lat- 
tice spacing. In order to estimate a systematic error in 
our combined chiral and continuum extrapolation, we 
choose fits, where the sea quark dependence is mod- 
elled according to the prediction of Heavy Meson Chi- 
ral Perturbation Theory (UMxPT) fSOl |5T1, as well as 
only linear in Wpg, 



/B(mps,fl,HYPn) 
b' 



(10) 



3 l+3g^ 2 , , 2 ^ 
1 - T -. ... '«ps In (wps) 



-H C 'Wpg -h d'i_ 



4 (47r/j2 

/b, (mps,a, HYPn) = b^ c^ mpg -H t/r,HYP« ; (1 1) 

here, f„ = f^^"^ and g = 0.51(2) are the same as above. 
These joint extrapolations are depicted as the black solid 
curves in Figs.[5]and[6] In particular for /b , not all CLS 
ensembles were analysed yet. 

From the figures one can infer that whether we do or 
do not include the chiral logarithm of HM;i^PT in the 
extrapolation ( [TOj l of /b induces a very small change at 
the physical point only. We thus take the HM;^'PT ex- 
trapolation as the central value and the difference to the 
linear fit to account for a part of the systematic error of 
our final result. For the B- and Bj-meson decay con- 
stants from HQET at 0(l/OTb) in two-flavour QCD we 
preliminarily give 



/b 

/Bs 



193(9)(4)^MeV, 
219(12) MeV, 



(12) 
(13) 



For the subtracted bare quark masses appearing here, the additive 
improvement term ii^^'amq j is numerically very small in practice. 
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HMxPT: /B(zb,mps,a)/GeV 



/[j(zb,mre,fl}/GeV 



/B,(zb, ;jire,fl)/GeV 




[mps/MeV]2 



Figure 5: Joint chiral and continuum extrapolation to the physi- 
cal point of the B-meson decay constant ^ in NLO HQET to the 
HMxPT-motivated ansatz \10\ . The colour coding is the same as in 
Fig-^ (I.e., blue, red and green points refer to p = 5.2, 5.3 and 5.5, 
while filled/open symbols belong to the HYPl/2 static actions.) 



where the quoted errors again cover all sources of sta- 
tistical and systematic uncertainties. 

Our results ([12]) and ( pj) are in line with computa- 
tions of other groups, see, e.g., summaries in ||52ll47l . 



4. Outlook 




Figure 6: Joint chiral and continuum extrapolation to the physical 
point of /b (left) and /b^ (right) in NLO HQET, where only a lin- 
ear dependence on the squared light pseudoscalar (sea) mass rrt^^ 
is assumed, cf. \1 In case of /bs . the NLO HQET result /b^ = 
216(5) MeV obtained in the quenched approximation (Nf = 0) fiSt 
(where the scale was set through rg = 0.5 im) is included for compar- 
ison. The colour coding is the same as in Fig.|5| 
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The non-perturbative treatment of NLO HQET with 
controlled chiral and continuum extrapolations leads to 
results for B-physics phenomenology with a few-% ac- 
curacy. They can contribute to resolving current ten- 
sions in precision CKM analyses of the B-meson sec- 
tor As our computations are the only ones, which have 
no perturbative uncertainties, including the renormal- 
ization of the axial current, this introduces a new quality, 
albeit for Nf-2. 

We are also investigating further spectral quantities 
within our approach, for instance, the B-meson spin 
splittings. Owing to the heavy quark spin-symmetry, the 
mass difference between the vector B*- and the pseu- 
doscalar B-meson is dominated by a pure 0(l/m[,) ef- 
fect from the contribution of (9spin to the effective HQET 
Lagrangian ([TJ and, hence, is of particular interest. 

By extending our finite-volume matching strategy to 
all components of the axial and vector currents, we aim 
at a NLO HQET calculation of B — » tt semi-leptonic de- 
cay form factors as possible application. A status report 
in the LO (static) approximation has been given in ||53]| . 
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